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ABSTRACT 

Young tidal dwarf galaxies (TDGs) are observed in the tidal debris of gas-rich 
interacting galaxies. In contrast to what is generally assumed to be the case for isolated 
dwarf galaxies, TDGs are not embedded in their own dark matter (DM) sub-halo. 
Hence, they are more sensitive to stellar feedback and could be disrupted on a short 
time-scale. Detailed numerical and observational studies demonstrate that isolated 
DM-dominated dwarf galaxies can have lifetimes of more than 10 Gyr. For TDGs 
that evolve in a tidal field with compressing accelerations equal to the gravitational 
acceleration within a DM sub-halo typical of an isolated dwarf galaxy, a similar survival 
time is expected. The tidal acceleration profile depends on the virial mass of the host 
galaxy and the distance between the TDG and its host. We analytically compare 
the tidal compression to the gravitational acceleration due to either cuspy or cored 
DM sub-halos of various virial masses. For example, the tidal field at a distance of 
100 kpc to a host halo of 10 13 M 0 can be as stabilizing as a 10 9 M 0 DM sub-halo. 
By linking the tidal field to the equivalent gravitational field of a DM sub-halo, we 
can use existing models of isolated dwarfs to estimate the survivability of TDGs. We 
show that part of the unexpectedly high dynamical masses inferred from observations 
of some TDGs can be explained by tidal compression and hence TDGs require to 
contain less unobservable matter to understand their rotation curves. 
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1 INTRODUCTION 

In the current standard cosmological theory, A-CDM, where 
the Universe is dominated by dark energy (A) and galaxies 
by their cold (C) dark matter (DM) content, dwarf galax¬ 
ies (DGs) are the remainders of the hierarchical structure 
formation process. They form the low-mass end of DM- 
dominated objects, as they have not yet merged into more 
massive systems. State-of-the-art cosmological simulations, 
such as Illustris (Vogelsberger et al.|2014|, EAGLE (Schaye 


galaxies in the Local Group (Holtzman, Afonso & Dol 


et al.|[2014| and zooms of Local Group analogues with the 


EAGLE code (Sawala et al.|2015 1 can now follow the evolu¬ 


tion of hundreds or thousands of DGs in one simulation box 
with increasingly high resolution. The stellar feedback pro¬ 
cesses that are especially important to regulate star forma¬ 
tion (SF) in low-mass galaxies, can be calibrated to observed 
relations, such as the galaxy stellar mass function. However 
not all discrepancies between simulations and observations 
on the scale of DGs are solved (Kroup a|2015 1. 

The simulations complement deep surveys 
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phin 

2006 Orban et al. 2008), the local volume (Hunter 

et al. 

2012 Lee et al.|j2008 

or in the Virgo (Lisker et al. 

2007| 

De Looze et al. {2013 

|Toloba et al. 2015} and For- 


nax clusters (Mieske et al. 2006 Chilingarian et al. 2011). 


Both observations and simulations indicate that very ex¬ 
tended (> 10 Gyr) lifetimes are possible for isolated DGs 
(e.g. Sawala et al.|2015 Karachentsev et al.|2014 2015 I. The 
explanation for their long-time survival is usually found in 
the DM sub-halo that engulfs the DG, increases the binding 
energy and stabilizing it against disrupting processes such 
as stellar feedback. 

Tidal dwarf galaxies (TDGs) form in the tidal debris of 
interacting gas-rich galaxies. In A-CDM, the merging host 
galaxies are surrounded by a DM halo, but the velocity dis¬ 
persion of these DM particles is too high to be captured by 
newly formed, low-mass TDGs which are therefore not ex- 



pected to be hosted by DM sub-halos l 

Barnes & Hernquist 

of dwarf 

1992 

Bournaud & Duc]|2006. Wetzstein, Naab & Burkert 


2007 

). Observations show that TDGs have active SF ( 

Braine 


the TDG on a short time-scale without a DM sub-halo to 
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keep it stable. Old TDGs are difficult to identify, as their 
prominent birth-place, the tidal arm, typically dissolves on 


a time-scale of several hundred Myr (e.g. Hibbard & Mi 


hos 1995 Bournaud & Due 20061. Without detailed kine¬ 


matic or spectroscopic studies, old TDGs might be hard 
to distinguish from other stellar objects in the same mass 
range. Only a few simulations of the evolution of individ- 


ual TDGs (e.g. Recchi et al.||2007 

Smith et al.|2013 Yang 

et al. 2014 Ploeckinger et al.|2014 

2015) or their formation 

during galaxy interactions (e.g. V 

Wetzstein, Naab & Burk- 

ert||2007; Bournaud, Due & Masset||2003; Bournaud & Due 

2006 Hammer et al.|2010 Fouquet et al. 2012) are available 


in the literature. 

While cosmological simulations are beginning to resolve 
the mass scales relevant for TDG formation, the feedback ef¬ 
ficiencies are sometimes taken to be functions of the proper¬ 
ties of the DM (e.g. Vogelsberger et al. 2014), which is not a 


good strategy for the stellar feedback inside TDGs. If many 
TDGs survive for several Gyr, then the number of old TDGs 
hiding in the DG and super-star cluster populations could 
be large. The first observational evidence for their long-term 
survival was found in a survey around early-type galaxies, 
including the detection of the oldest TDG identified so far, 
with an age of around 4 Gyr (Due et al.|2014|. 


In Ploeckinger et al. (2015) the evolution of individ¬ 


ual TDGs is simulated for 3 Gyr. During their evolution, 
the TDGs pass both the apo- and the peri-centres of their 
eccentric orbits around their host galaxy. While they are 
close to the apo-centre, the star formation self-regulates to 
a moderate level, but the TDGs become significantly tidally 
compressed as they get near the host galaxy]^] Unlike star 
clusters, where tidal compression leads to dynamical heating 
of the stellar system, TDGs are gas-rich and dissipative pro¬ 


cesses are important for their formation (Wetzstein, Naab 
|fe Bur kert 2007j). The tidal field compresses the gas phase, 
leading to efficient cooling and a high star formation rate. 
In addition, this compression helps to stabilize the TDG 
against the subsequent stellar feedback (fPlocckingcr et al. 


2015). Kroupa (19971 and Klessen & Kroupa (19981 have 


demonstrated that gas-free TDGs can survive for 10 Gyr 
despite repeated tidal harassing. 

Tidal forces are often used to study the disruption of 
objects (such as stars close to a black hole), but a tidal field 
can also trigger the formation of new systems or increase 
their stability in its compressive mode (jRenaud et al.|2009| 


2014). In regimes where the compression of the tidal field 


dominates, an object experiences an extra acceleration to¬ 
wards its centre, reminiscent of an additional gravitational 
acceleration as e.g. caused by a DM sub-halo. By comparing 
the acceleration due to the tidal field to that caused by a 
DM sub-halo of a given mass, we can relate the evolution 
of TDGs to that of their much better studied isolated DG 
siblings: If an isolated DM-dominated DG of a given mass 
can survive its stellar feedback for almost a Hubble time, 
then a similar TDG at a position within a tidal field with 
comparable accelerations should survive as well. 

We present a fully analytic method to assign an equiva- 


1 For a movie of the tidal compression of model TDG-p from 
[Plo eckinger et ah] ( j2015[ l , see https://sites.google.com/site/ 
sylviaploeckinger/tidal-dwarf-galaxies/analytical 


lent DM sub-halo mass to a TDG depending on the accelera¬ 
tion due to the tidal field and therefore on the mass and pro¬ 
file of their host halo and the distance between the TDG and 
the centre of the host galaxy. The method provides a simple 
tool to predict or explain the survival of TDGs. Bournaud 


et al. (2007) measured the rotation curves of three TDGs in 


NGC 5291 and found the ratio between the dynamical and 
the visible mass to have a value of 2 to 3. As TDGs should 
be free of DM, these ratios are unexpectedly high. We show 
that the compression of the tidal field can lead to a com¬ 
parable increase of the dynamical mass as a DM-halo with 
a mass equal to the “missing mass” in these TDGs could. 
Therefore this method can help to explain an apparent DM 
content in TDGs. 

In Sec. [2] the general properties of tidal fields are sum¬ 
marized and applied to find distance limits for each host halo 
mass, in which the tidal compression along two directions 
dominates the tidal stretching in the perpendicular direc¬ 
tion. In Sec. [3] the accelerations from the tidal compression 
are compared to those inside DM sub-halos. In Sec. [4] we 
present two applications of this method: explaining the age 
of the old TDG in NGC 5557 and investigating the high 
dynamical masses of the TDGs in the gaseous ring around 
NGC 5291. We discuss the limitations of the proposed ana¬ 
lytic method in Sec. [5] 


2 TIDAL FIELDS 

For the method presented here, we consider a host galaxy 
centred at x = y = z = 0 with a spherically symmetric 
gravitational potential 4>host from a mass distribution phost 
and investigate the acceleration field at the surface of a test 
sphere with radius r at a distance D, as shown in Fig. |XJ 
The acceleration field in the reference frame of r = 0 
(tidal field), caused by the external potential is described as: 


= V$ hos t(0) - V$ host (r') 


(1) 


where r' is the position vector relative to the centre of the 
test sphere. The gravitational potential of the host galaxy 


is described by an NFW profile (Navarro, Frenk & White 


19951 of 


Pnfw(R) 


R 

r 0,h 


pO,h 




with the resulting potential 


( 2 ) 


4?nfw (R) = — 


Gpo,h.r oh 

R 


log(l + R/r 0 , h ) 


(3) 


where po,/i is the characteristic density of the host halo and 
ro : h its scale radius. The tidal field at z = 0 for a test sphere 
of r = 3 kpc at a distance of D = 10 kpc to the centre 
of an NFW halo with M v i r = 10 12 Mq is indicated as vec¬ 
tors in Fig. [l] for illustration. In this case, the test sphere 
would be stretched in ^-direction and compressed along the 
y- and z axis. In a potential with a central cusp, such as the 
NFW profile, the tidal field is never compressive on all axes 
(see e.g. Renaud et al.~| |2009 j. Nevertheless, it is possible to 
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Figure 1. The vectors show the tidal field at a distance of 
D = lOkpc to the centre of an NFW potential at x = y = 0 . 
Throughout the paper the definitions for x, y, D, r = |r'|, R, o . r , 
and a y are used as indicated here. 



0 50 100 150 200 


D [kpc] 

Figure 2. The colour scale represents the ratio between the eigen¬ 
values of the tidal tensor A 2 /i at distances D (x-axis) to the cen¬ 
tre of NFW potentials with scale radius ro (y-axis). The solid 
black line shows where A 2 /i = —1 (see Eq. |10l. The blueish re¬ 
gion above the solid line corresponds to a contraction dominated 
regime (A 2/ q < —1), while in the reddish region below the black 
line, the expansion dominates (A 2 /j > —1). For more massive ha¬ 
los and therefore NFW potentials with larger scale radii r o, the 
region in which the compression dominates the expansion extends 
out to larger distances. For more information see text. 


look for regimes, where the compression dominates the ex¬ 
pansion, which is the aim of this section. Note that during 
a galaxy interaction, the gravitational potential will not be 
relaxed. In principle, the method can be used for arbitrary 
gravitational potentials if the analytical are replaced by the 
numerical derivations. For simplicity, we show the idealised 
case of an NFW potential, but see the discussion on this in 
Sec. [5] 

2.1 Gauss’s law of gravity 

The net flow over the whole surface of the test sphere, shown 
in Fig.[l]is directly given by the integral form of Gauss’s law 
for gravity: 


host galaxy leads to an increasing enclosed DM mass and 
therefore to a higher net compression. 

2.2 Tidal tensor 

One way to describe the tidal field is with the tidal tensor 
TT = —di(dj‘ Fhost)- TT is a symmetric tensor and in its or¬ 
thogonal form, the eigenvalues A represent the strength and 
direction of the tidal field along the corresponding eigenvec¬ 
tors v. If A i < 0, the tidal acceleration leads to a compression 
in the direction of Vi while Ai > 0 represents an expansion 
along Vi. 

In the geometry illustrated in Fig. [I] the components 
of the tidal tensor at a distance D to the host galaxy are 
defined as: 


g ■ dA = —AnGMe, 


(4) 


It states, that the integral of all accelerations g in a gravi¬ 
tational potential over any static surface of the volume dV 
is proportional to its enclosed mass M enc . For the above ex¬ 
ample this means, that if the massless test sphere is outside 
of the density distribution phost and therefore M enc = 0, the 
net effect over the surface elements A vanishes. The tidal 
stretching in some directions is therefore exactly balanced 
by the tidal compression in other directions. If the satel¬ 
lite object is inside the mass distribution, then M enc > 0 
and subsequently the net effect on a static surface is always 
compressing [§ gv g ■ dA < 0 ). 

Therefore, whenever a TDG is inside a DM host halo, 
there is a net compression caused by the tidal field. In a 
related work, Bournaud, Due & Masset (20031 showed that 
massive TDGs at the tips of tidal arms only form in their 
simulations if the DM halo of the host is not truncated. 
Outside the truncation radius, the tidal field doesn’t lead to 
the necessary compression for the formation of TDGs. 

A first condition for the formation and survival of TDGs 
is therefore to reside within the host galaxy’s DM halo which 
translates to D < r v n, where r v ir is the virial radius. In 
addition, the net compression over a surface is proportional 
to the enclosed mass within it. A decreasing distance to the 


T 11 (D) 

T 22 (D) 

T 33 (D) 


/cohost (D)\ 

v dx 2 ) 

/ cohost (D)\ 

V dV 2 J 

fd 2 $ host (D)\ 

V dz 2 ) 

0 


(5) 

( 6 ) 

(7) 

( 8 ) 


with the eigenvalues Ai = T 11 , A 2 = T 22 , and A 3 = T 33 ) and 
the eigenvectors Vi = ( 1 , 0 , 0 ), vn = ( 0 , 1 , 0 ), and V 3 =( 0 , 0 , 1 ). 

The trace of the tidal tensor can then be directly linked 
to Poisson’s law: 


3 3 

A* = I] - 8 i(ft*hart) = -V 2 $ho S t = —4 t rGp . (9) 

i=l i=l 

As the density p is always positive there is no solution where 
all A i are positive and where therefore the tidal field would 
lead to an expansion along all eigenvectors. Depending on 
d'host and the distance D, between one and three eigenvalues 
can be negative and therefore compressing the test sphere in 
the corresponding direction. In the example shown in Fig. [l] 
the eigenvalue in x-direction is positive while those in y- 
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D [kpc] 


ro 

[kpc] 

M v ir 

log [M/Mg] 
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Figure 3. The ratio between the tidal acceleration in y- and x-direction {a y / x , see also Fig. |Tjl is presented in dependence on the distance 
to the host galaxy (x-axis),the radial extent of the TDG (y-axis), and the scale radius ro of the host halo (solid lines). All lines show 
a y/x = 1 f° r NFW potentials with different scale radii according to Eq. |12| and are labelled with the value for ro in kpc. The grey regions 
indicate where a y / x > 1 and therefore the compression dominated region for each scale radius. The corresponding virial mass M v j r for 
a halo with scale radius rg is listed in the table on the right (see Sec. |2.3| lor details). 


and z-direction are negative. Fully compressive tidal modes 
where all Ai < 0 appear in central regions of cored poten¬ 
tials or in the overlay of two potential e.g. as in interacting 
galaxies. For a detailed description of the tidal tensor and a 
discussion on the different modes of a tidal field, see|Renaud| 


et al. (20091. 


Even if the tidal mode is not fully compressive, it’s effect 
can still be more stabilizing than disrupting. An indication 
for that is the ratio between Ai and A 2 . The fully compressive 
mode is then described when all Ai < 0 while for A 2 /A 1 < — 1 
the contraction in y-direction is stronger than the expansion 
in x-direction. Note that for the geometry used here, the 
problem is symmetric around the x-axis and therefore A 3 = 


A 2 . 

For an NFW potential (Eq. |3j) the ratio A 2 /i = A 2 /A 1 
can be calculated analytically]^] and depends only on D and 
the scale radius r q: 


A 2 /i (D,r 0 ) = 

Z? 2 + ,Dro-(D + ro) 2 log(T^o) (10) 

2 (D 2 + r 0 (2 D + r 0 )) log - 3 D 2 - 2 Dr 0 

In Fig. [2|the relation between D, ro and A 2/1 is illustrated. 
The black solid line indicates where the contraction in y- 
(and z-) direction is equal to the expansion in x-direction. 
Above the A 2/1 = —1 line, the contraction is larger than the 
expansion, leading to a first limit on where the survivability 
of objects could be enhanced rather than decreased. 


potential of the host, but also on its radial extent r. As it 
can be directly seen from Eq. [T] the more extended the test 
sphere is, the larger is the difference between the acceleration 
at its mass centre and those at its surface. 

In addition to the dependence of the absolute values of 
the tidal forces on the radius r, also the ratio between com¬ 
pression and expansion changes with r. In order to include 
the radial dependence, we calculate the ratio a y / x = \a y /a x \. 
We define a y as the acceleration along the y-axis at the point 
{x = D,y = r,z = 0; see Fig. [I]) and use it as a proxy for the 
maximum compression. Analogously, a x presents the accel¬ 
eration in x-direction at the point [x = D — r, y = 0, z = 0) 
and is the proxy for the maximum expansion. 

Compared to A 2/1 that depends only on D and ro, a y / x 
also provides information on the different modes within the 
test sphere or for test spheres with different radial extents. 

The y-component of the tidal field at the position ( D , r) 
in the reference frame of the host galaxy (see a y in Fig. [l| 
is given by 2 : 


Kl 



d$(R) 

OR 


(ID 


where R = yjD 2 + r 2 is the distance of ( D,r ) to the host 
galaxy. For an NFW potential, a y / x is therefore: 


2.3 Tidal accelerations 

The tidal field at the surface of the test sphere does not only 
depend on the distance to and the shape of the gravitational 

2 The Mathematica notebook for the derivation and further 
use of Eqs. |10| |12| and |22| as well as an illustration for 
Eq. |lf| can be downloaded at https://sites.google.com/site/ 
sylviaploeckinger/tidal-dwarf-galaxies/analytical 


a y / x (r,D,r 0 ) = 

(Vr 2 + D 2 + r 0 ) log ^ \A 2 +^ 2+r o ^ _ A / r 2 + j >2 ( 12 ) 

(■ r 2 + D 2 ) sj^r + l (V r 2 + D 2 + r 0 ) {A + B) 

with: 
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(~r + D + r 0 ) log ( - r+ ” +r ° ) + r - D 
(r — D) 2 (—r + D + r 0 ) 

D-(D + r 0 )lo g (°^) 

D 2 (D + r o) 


(13) 

(14) 


For a scale radius ro of the external NFW potential, the 
regimes where the compression a y is larger than the expan¬ 
sion a x are defined by a y / x > 1. Fig. [3] illustrates these 
regions for scale radii ro between 10 and 100 kpc. For larger 
values of r o and therefore more massive host galaxies, the 
regions where compression dominates get bigger, they ex¬ 
tend out to both larger distances D as well as to larger radii 
r. 

The virial mass M V - 1T for a DM halo is defined as 


with 


and 


A/vir = 47rp 0 ro ( log(l + c) - 


pO — per t) c 


1 + C 


(15) 


(16) 


Sc = 


200 


(17) 


3 log(l + c) — c/(l + c) 

where p CI = 3Ho/(8nG) is the critical density of the Uni¬ 
verse, while Sc is the characteristic over-density of the in¬ 
dividual DM halo and dependent on the concentration pa¬ 
rameter c. Correa et al. (2015) derive the following fitting 
functions for the c—M v i r relation (see Sec. [5] for a discussion) 
in the Planck cosmology: 


2011; Amorisco &; Evans||2012 1 . Baryonic processes such as 


stellar feedback and a resulting wind are proposed as a mech¬ 
anism to produce DM profiles with a central core (see e.g. 


Pontzen & Governato 20121 but its realization is not yet 
fully understood ( Recchi &; Hensler|2013 1. In this section we 
compare the tidal field inside a massive NFW halo to DM 
sub-halos with both cuspy (NFW-NFW) and cored central 
regions (NFW-Burkert). As an example for a cored DM dis¬ 
tribution, we use the Burkert halo, an empirical profile from 


Burkert (19951: 


PBurkert (A) — 


Po»"o 


(R + ro)(R 2 + rg) 


(19) 


with the central density po and a scale radius ro. 

The accelerations within each DM profile are shown in 
the first two panels of Fig. [ 4 ] for sub-halo masses of 10 s , 
10 9 , 10 1 0 , and 10 11 Mq. The third to fifth panel of Fig. [i| 
demonstrate the accelerations along the y-axis (see Fig. fTF 
for host masses of 10 11 , 10 12 , and 10 13 ,Mq for distances of 
30, 60, and 90 kpc. 

Similar to a y / x , we can compare the accelerations 
at, y (y) to aoM(r), and find regimes where the tidal acceler¬ 
ation dominates and therefore a t /DM = at, y {y)/ar>M(r) > 1 

For each distance D to the host galaxy and for each 
radius r of a test sphere in this tidal field, the compressing 
tidal force can be related to the acceleration within a DM 


sub-halo with either a central cups (NFW, Sec. 3.1 1 or core 
(Burkert, Sec. |3.2|>. 


3.1 NFW tidal field compared to an NFW 
sub-halo 

We are comparing here the compressing tidal accelerations 
of an NFW host halo a y to the gravitational accelerations 
inside an NFW sub-halo o r : 


logio c = a + /?log lo (Af/M 0 ) [1+ 7 log 2 o (Af/M 0 )] (18) 

with a = 1.49809, /3 = -0.02499, and 7 = 0.0056 for z = 0. 
As pc and c depend on A/ V i r in a non-trivial way, we 



3 COMPARISON TO DARK MATTER 

SUB-HALOS 

In the previous section we have shown at which positions a 
tidal field can lead to a compression in y- and z-direction 
that is larger than the expansion it causes in x-direction. 
In its astrophysical application, this additional force can 
enhance star formation but could also stabilize the object 
against stellar feedback, similar to a DM sub-halo. 

Here we compare the accelerations inside a DM sub¬ 
halo a r = \d<&/dr\ to the compressing acceleration caused 
by the tidal field a y (Eq. ED- Pure N-body simulations pre¬ 
dict NFW profiles for DM halos on all scales, but obser¬ 
vations show indications for centrally cored density profiles 
for DCs (e.g. |Kuzio de Naray, McGaugh fe de Blok 2008 


|Battaglia et ah] 2008| |De Blok| 2010| |Walker fe Penarrnbia| 


47rGpo, t ro |t ( r Q4 + R ) lQ g ( r 0|t ) R 

^+(f) 2 


R 2 (r 0l t + R) 


( 20 ) 


where R = V D 2 + r 2 as in the case of a y / x (Eq. 12 ) and 


= ATvGpo,d.ro, d 


(r 0 ,rf + r) log (~~~) 
r 2 (r 0 ,d + r) 


( 21 ) 


pa t t and ro,t describe the potential of the host galaxy for the 
tidal field while po.d and ro,d are the characteristic density 
and the scale radius of the DM sub-halo. The ratio a y / r = 
|a^/a r | between the tidal compression and a DM sub-halo 
can be calculated with 2 : 


where 
A = 

B = 


Po,t ro t 1 A 

ay/r=p °’ d 7(|FtT b 

r 2 (r -(- r 0 ,d) R + r 0 ,t) log 
R 2 {R + ro,t ) ^(r + r 0 ,d) log 


( 22 ) 


R) (23) 
r) (24) 
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NFW Burkert Tidal field Tidal field Tidal field 



r DM [kpc] r DM [kpc] r ty [kpc] r ty [kpc] r ty [kpc] 

Figure 4. Comparison of the accelerations inside a DM sub-halo with an NFW (first panel) and a Burkert profile (second panel) as well 
as those in a tidal field at distances 30 (third panel), 60 (fourth panel), and 90 kpc (fifth panel) to the host galaxy. The accelerations are 
indicated for different sub-halo virial masses M su b (first two panels) and host masses A/host (last three panels, see legend). We assume 
an NFW potential for the host galaxy in the case of the tidal field. Unlike the DM potential, the tidal field is not spherically symmetric 
around r = 0. We use here the compressing acceleration a y as indicated in Fig. |TJ along the y-axis. In order to highlight the asymmetry 
the radius for the tidal field is labelled as rt lV . 


Assuming the relation between the concentration pa¬ 


rameter c and the virial mass (Eq. 181 from Correa et al. 


120151 allows to fully describe each halo profile with only 
one parameter, here the virial mass M v n. The characteristic 
density po follows then from Eq. [16] an d|17[ while the scale 
radius can be calculated by solving Eq. 1 15 1 for ro- 

The relation a y / r therefore depends on r, D, and the 
virial masses of both the host galaxy and the comparison 
DM sub-halo. Fig. [5] shows the lines for a y / r = 1 for various 
host and sub-halo masses. The lines in each panel represent 
the edge of the regime where the compression of the tidal 
field is stronger than the DM sub-halo it is compared to. As 
the tidal field increases with larger r while the acceleration 
inside an NFW sub-halo decreases at larger radii (see Fig. [4j) , 
the area above each line (for r > r[a y / r = 1)) represents the 
regime where the tidal field dominates. 

An example: the compressing tidal field at a distance 
of D = 100 kpc to the center of a host halo of M v i ri ho B t = 
1O 13 M 0 (Fig. bottom panel), is equivalent to the acceler¬ 
ation of a DM sub-halo of M v i r , S ub = 10 8 Mq at a radius of 
r = 2.5 kpc. At larger radii, the tidal field gets stronger and 
it therefore compares to a DM sub-halo of A/ v j rjSU b = 10 9 M 0 
at a radius of 5.5 kpc and M v i r , su b = 10 10 M 0 at r = 11 kpc. 
For applications and interpretations see Sec. [4] 


3.2 NFW tidal field compared to a Burkert 
sub-halo 

Analogously to Sec. |3.1| we repeat the same analysis for DM 
sub-halos with a Burkert profile. As we assume an NFW 
profile for the massive host galaxy, a y is still represented by 
Eq.|TTj The radial acceleration inside a Burkert halo is given 
by: 


7/3 

7V(jr Pconst^Q 


2 arctan —J + 41og(r 0 ) — 21og(r + r 0 ) - log(r +ro) 

(25) 

where po is the central DM density, p C onst a constant that 


includes all constant factors for po in cgs units (see Eq. 271, 
and ro is the scale radius of the Burkert halo (Eq. |19[). Ac¬ 


cording to Burkert (19951, they correlate with each other 


and the virial mass M v i r via: 


po = 4.5 x 10“ 


r o 


kpc 


-2/3 


M 0 pc 


-2/3 -3 

Pconst r 0 gem 


Mv 


= 5.8 x 7.2 x 10' x 


r o 
kpc 


-i 7/3 


M 0 


The results for a y / T = 1 are shown in Fig. [6] 


(26) 

(27) 

(28) 


4 APPLICATIONS 


Cosmological simulations are now able to resolve dwarf 
galaxies with unprecedented details where effects of the vari¬ 
ations of stellar feedback processes are studied carefully by 
different groups, as e.g. the EAGLE ( Schaye et al.|2014 ), the 
FIRE ( Hopkins et al.|2014 1, or the Illustris project ( |Vogels- 
berger et al.||2014 |. All of these objects are DM-dominated 
within A-CDM. 

The data in both simulations and observations on young 
TDGs is very limited compared to the DM dwarf galaxies. 
Here, we relate the properties of the tidal field around TDGs 
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Figure 5. Comparison between the accelerations inside NFW DM sub-halos with virial masses M v j r gut, of 10 s (red solid lines), 10 9 
(blue dashed lines), !0 19 (green dotted lines), and 10 11 M... (grey dash-dotted lines) and in a tidal field of a host galaxy with virial 
masses M v ; r host from 10 11 to I (1 1,1 M... (top to bottom panel). All lines represent the border of the region above which the tidal field is 
stronger than the representative DM sub-halo, therefore where Ut/DM = 1 (Eq. |22| . For larger radii r or smaller distances to the host 
D, the tidal field dominates. The grey shaded area represents the region where the compressing acceleration in y direction is larger than 
the expansion in x-direction (a y ^ x > 1. compare Fig.|3l. 


to those of the DM halo of isolated dwarfs. This way it is 
possible to estimate the fate of TDGs as well as explain 
the peculiarity of the apparent DM content in some of the 
TDGs. 

This exercise assumes perfect NFW potentials for both 
the host galaxy as well as the DM sub-halo. This has to be 
taken with caution, especially as the potential of the host is 
time-variant during the galaxy interaction. Nevertheless it 
provides a quick and easy check on the acceleration regime 
in which the TDG is located, and allows first estimates on 
the potential survivability of TDGs. 


4.1 NGC 5557 


Due et al. (20141 found three TDGs around the early-type 


galaxy NGC5557. SED fitting of the TDG NGC5557-E1 
reveals an age of 4 Gyr, which makes it the oldest TDG 
uniquely identified so far. The work of this paper provides 
a possible explanation for the long-term survival of this ob¬ 
ject._ 


Cappellari et al. (2013) derived a maximum value of the 


circular velocity for NGC5557 of v max = 340kms _1 from 
their best-fitting mass-follow-light JAM (Jeans Anisotropic 
Multi-gaussian expansion) model. The virial mass of NGC 
5557 can therefore be estimated with 


M vir = 2.33 x 10 


5 / ^max 

Vkms~ 


0‘ 


Mp, 


(29) 


for an NFW potential (Navarro 1998). For the TDG 


NGC5557-E1 the host mass is therefore in the order of 
log 10 A/ Vj ngc5557 = 13. This is a very crude mass estimate, 
but to illustrate the application of this method, an order of 
magnitude approximation is enough. 


From the images in Due et al. (2014), the projected 
distance of the TDG to its host is about D as 70 kpc in 
projection and the HI contours show a spatial extent be¬ 
tween r « 5kpc (minor axis) and r as 6kpc (major axis). 
The combination of log M v i ri host [Mq] as 13, D > 70 kpc, and 
r = 0 — 6 kpc can be overlaid to Fig. [5] which is shown as a 
short black vertical line in Fig. [7] 

At the edge of TDG NGC 5557-El (r as 6 kpc Due et al. 


2014) the accelerations are comparable to those of an NFW 
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Figure 6. As Fig. [ 5 ] but here the DM sub-halos are represented with a cored potential from|Burkert| (1 995[- 


DM sub-halo with a few times 10 9 Mq or to a Burkert sub- 
halo with 10 9 Mq. Around the effective radius of NGC 5557- 
El ( r e = 2.3kpc), where the baryonic matter is expected to 
dominate in all cases (Toloba et al. 2011), the additional 
acceleration from the tidal field is still the same as in DM 
sub-halo with M v i rjSU b ~ 10 s Mq. In addition, TDG NGC 
5557-El is also well in the regime where the compression in y 
and z direction dominates the tidal stretching (grey shaded 
area in Fig. [7]). 

This might be an explanation for the extended lifetime 
of this particular TDG. Detailed kinematical studies of NGC 
5557 and its surroundings, used for a DM halo profile fit, 
would allow for a much better estimate. In principle, once 
the gravitational potential of the host galaxy is known, these 
regions can be constructed easily. 

As an estimate for the further evolution of this object, 
the tidal radius can be estimated with 


_ f Mtdg \ 1/<3 
' j ~ \3M host (D)J 


(30) 


(Binney & Tremai ne[20 08), where Mhost(D) is the enclosed 
mass of the host halo up to the distance D. For NGC 5557- 
El with a mass of Mtdg = 1-2 x 10 8 Mq and at a distance 
of 70 kpc, the r, is at 2 kpc, close to the current r e . 

As the estimates on the mass and DM profile of NGC 


5557 are only based on one measurement of the maximum 
rotation velocity, no precise statements can be made yet 
about the individual objects. But some general conclusions 
for a TDG at a distance of 70 kpc to a host galaxy with a 
virial mass of 10 13 Mq and an NFW profile can be drawn 
nevertheless: 


- Around 5 kpc away from the center of the TDG, gas 
from the surrounding tidal arm is compressed by the tidal 
field onto the TDG, as if the TDG was embedded in a DM 
halo of a few times 10 9 Mg. 

This effect increases with increasing radii (10 1o Mq at 
11 kpc for the NFW-NFW case), which is an efficient way 
to bring gaseous material close to the centre of the TDG. 

- The tidal field can fuel SF and stabilize the outer parts 
of the TDG against stellar feedback, just as a DM halo in 
this mass range could. 

I 11 the very central region of the TDG, the contribution 
of the tidal field vanishes. But also DM-dominated DGs are 
dominated by the self-gravity of the baryonic matter close 
to their centre. 

- At this distance (70 kpc), the TDGs are clearly in the 
regime where the compression of the tidal field dominates 
the tidal stretching (a y / x > 1). 

For a TDG with 10 8 Mq the dynamics at the outer parts 
could mimic a DM content in the order of a few times the 
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Figure 7. Labels as in Fig. [5] In addition the estimated distances of the TDGs around NGC 5557 and NGC 5291 are indicated as black 
vertical lines at D = 70 kpc (NGC 5557-El) and D = 130 kpc. For the NGC 5557 system, D is the projected distance and therefore 
a lower limit, as indicated by an arrow. |Bournaud et ah] ( |2007| performed simulations of the galaxy interaction that caused the tidal 
features around NGC 5291 and their best model places the TDGs at a distance of around 130 kpc today. Because of the special geometry 
of the interaction, they spent most of their evolution closer to the host galaxy than it is observed today. This is indicated by an arrow 
towards smaller D. The top panel shows the case NFW-NFW (Sec. E3 and the bottom panel the case NFW-Burkert (Sec. |3.2[). 


baryonic mass, caused by the additional acceleration of the 
tidal field. 

- An estimate of the tidal radius gives r, = 2 kpc. Even 
if all the gas is blown out eventually, a bound stellar object 
with a radius of around 2 kpc could survive, if it’s on a 
circular orbit. 


4.2 NGC 5291 


Bournaud et al. (2007) measured the rotation curves of three 


TDGs in the gaseous ring of NGC 5291. They found dy¬ 
namical masses that are between two and three times larger 
than their visible masses. TDGs cannot contain bound DM 
particles as they originate from disk material and they are 
not massive enough to capture the high velocity dispersion 
particles from their host galaxies. If the discrepancy be¬ 
tween the visible and dynamical mass is real and not caused 
by e.g. projection effects, it has to be explained otherwise. 
Among the proposed solutions are a modified Newtonian 
gravity (|Ge ntile et al. |[2007| )), a molecular gas content that 
is 10 times larger than expected with the traditional CO to 
molecular gas conversion factor ( Bournaud et al.|2007 1, or a 
velocity anisotropy caused by the tidal field (|Kroupa|1997|). 

NGC 5291 is included in the HIPASS (HI Parkes All Sky 


Survey) database and Koribalski et al. (2004) measured a 
50% velocity line width of w$o = 637 km s -i . Assuming that 
the maximum rotation speed Umax is close to 0.5 x W 50 , the 
virial mass of NGC 5291 is estimated to be M V} NGC5291 = 
7.5 x 10 12 M 0 (Eq. 291 or log 10 M v ,ngc529i = 12.9. 

Simulations of galaxy interactions to create the mor¬ 
phology of the gaseous ring around NGC 5291 favour a sce¬ 
nario of a high speed collision between NGC 5291 and the 
massive elliptical galaxy IC 4329 around 360 Myr ago (Bour 
naud et al.||2007 ). Their best model has the ring located in 
a plane, 130 kpc away from the host galaxy today. 

For an estimate on the future of the TDGs around NGC 
5291, they are indicated in Fig.[7]at D = 130kpc. Their ro¬ 


tation curves were measured out to 5 kpc (Bournaud et al. 


2007), where their acceleration by the tidal field is compa¬ 
rable to that of a DM sub-halo with around 10 8 M 0 . This 
could account for some of the missing mass for the TDGs in 
NGC 5291. 


As the ring was produced by an high speed encounter, 
a significant part of the evolution of the TDG was closer to 
the host galaxy (indicated with an arrow in Fig. [tJ) . Never¬ 
theless, the TDGs today are on the edge of the regime where 
a y /x < 1 and therefore where the expansion dominates the 
contraction. If they are not on eccentric orbits that will bring 
them closer to the host galaxy, they are more likely to get 
disrupted. 

As in the case of NGC 5557, only little is known about 
the detailed DM profile of NGC 5291. We assume again an 
NFW potential, but if the disturbance that caused the tidal 
ring has happened only a few hundred Myr ago, the DM halo 
of NGC 5291 is very likely not re-virialised. In addition, the 
estimates of M v> NGC5291 and the distances are first order 
approximations and should be used with great care. More 
detailed data on the DM distribution from simulations or 
observations would allow stronger statements on their future 
survivability. 

In general, the following information can be extracted 
for a TDG at a distance of 130 kpc to a 10 13 M 0 host galaxy 
with an NFW DM halo profile: 


-At D = 130 kpc, the tidal field is much weaker than in 
the previous example wit h D = 70 kpc. 

The tidal radius (Eq. 301 for a an object with 5 x 10 8 M 0 
is r.j = 5 kpc. 

The TDG is located at the edge of the a y / x > 1 region. 
At a radius of 5 kpc, the tidal stretching in x-direction is 
even stronger than the compression in y- and z- direction. 
In this regime, the formation of TDGs is not enhanced by 
the tidal field, as less mass can be channelled from the tidal 
arm to the close vicinity of the TDG. 
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The compressing accelerations from the tidal field are 
comparable to those of a 10 s M@ DM halo, therfore the dy¬ 
namical stage of the TDG could mimic a low-mass « 10 8 M@ 
DM halo. 

If the TDG moves to larger distances, the tidal radius 
would increase while the compressing tidal field decreases. 
An object like that would be unsupported by any stabilizing 
extra acceleration field and if it is still gas-rich at this stage, 
stellar feedback may be efficient in expelling the gas. 

If the gravitational potential of NGC 5291 has a higher 
concentration parameter c than the average value of |Cor-| 
|rea et al. |2015| the tidal field is stronger, resulting in a 
higher dynamical mass in its TDGs. In order to explain an 
additional dynamical mass of ~ 2 x 10 9 Mq, as measured 
in |Bournaud et aTT| ( |2007| ) for the TDGs around NGC 5291, 
the concentration parameter of the DM halo of NGC 5291 
would have to be unrealistically large (^> 1000). 

Note, that rather than making clear predictions on indi¬ 
vidual TDGs, those objects were chosen to illustrate which 
information can be extracted with this simple approach. 


5 DISCUSSION 


We have presented a simple analytic method to investigate 
the survivability of TDGs by comparing the accelerations 
caused by the tidal field to those inside a DM sub-halo. 
In this section we discuss the limitations of the presented 
approach. 

In general, the compression due to the tidal field in¬ 
creases with decreasing distance to the host galaxy (Eq. |4j). 
Using our terminology, this translates into higher equivalent 
DM sub-halo masses and therefore an increased stability of 
the TDG against stellar feedback processes. The lower limit 
on the distance between TDG and host galaxy is considered 
by a calculation of the tidal radius. While the compression 
increases with decreasing distance, the volume within which 
test particles are gravitationally bound to the satellite ob¬ 
ject, rather than to the host galaxy, decreases. This volume 
is characterised by the tidal (or Jacobi) radius (Eq. |30[ ). We 
do not directly include a lower distance limit in the method 
here, but we compare the tidal radius to the physical extent 
of the object for the presented applications as an indepen¬ 
dent test. 

Another process that limits the survivability of TDGs 
with smaller distances D to the host galaxy is the orbital 
decay caused by dynamical friction ( Chandrasekhar||1943 1. 
Simulations of interacting galaxies have shown that TDGs 
that form close to the host galaxy (< 50kpc) can be re¬ 


accreted on a time-scale of a few hundred Myr (e.g. Hibbard 
fe Mihos|1995 Bournaud &; Duc|2006 l. Note, however, that 
the dynamical friction time-scale for a TDG is much longer 
than for a DM-dominated satellite galaxy with the same 
baryonic mass, as diTorb/d t is proportional to the total mass 
of the satellite ( Chandrasckh ar|l943[ ) . 

Summarizing, we expect a “sweet spot” at particular 
distances to the host galaxy, where the survivability of TDGs 
is enhanced by the tidal field. The distance range depends 
on the mass and shape of the DM halo of the host and 
the analytic method presented here can help to determine 
this range. It also allows one to predict the survivability of 
individual TDGs given their distance to their host galaxies 


and the masses of their host halos, if the survivability of 
DGs in DM halos with equal accelerations is known. The 
long-term stability of a TDG is favoured if the distance to 
the host galaxy is sufficiently small for the tidal field to 
stabilize it against stellar feedback, but also large enough 
that its tidal radius does not get too small and dynamical 
friction will not lead to orbital decay on a short time-scale. 

We made several assumptions and simplifications: 

Eccentric orbits: The method presented here only de¬ 
pends on the current distance. To predict the further evo¬ 
lution or to reconstruct the past, the orbit and therefore 
the time-dependence of D has to be taken into account. A 
highly eccentric trajectory will bring the TDG through very 
different regimes. Note however that long-living TDGs are 
expected to move preferentially on orbits with mild eccen¬ 
tricity as otherwise they would get too close to the host 
galaxy eventually. 

Host halo potential profile: We assume spherically sym¬ 
metric NFW halos for all host galaxies. At the late stages of a 
galaxy interaction, especially after the galaxies have merged, 
this assumption is expected to be valid. During the galaxy 
interaction event, the DM distribution is not virialised and 
might therefore deviate significantly from the NFW profile. 
Simulations can provide time-dependent gravitational po¬ 
tentials and the tidal compression can then again be com¬ 
pared to accelerations in DM sub-halos. 

c-M relation: For all NFW halos, including both the 
host galaxy and the DG, we assume the relation between 
the concentration parameter and the virial mass from |Cor^| 


rea et al. (2015). Their relation combines an analytic model 


for the halo mass accretion history based on the extended 
Press Schechter theory with a fit between the concentration 
parameter and the formation time from DM-only simula¬ 
tions. We consider the DM profiles of average host galaxies 
but real galaxies are expected to scatter around the c-M 
relation. 

Asymmetry of the tidal field: In contrast to the gravita¬ 
tional field inside an isolated DM sub-halo, the tidal field is 
not spherically symmetric. When we compare the accelera¬ 
tions, we take the maximum compression of the tidal held 
a y (see Fig. [lj. The extension in the direction towards the 
centre of the host halo is accounted for when the ratio be¬ 
tween compression and expansion is calculated. In the cho¬ 
sen geometry of this toy model, the tidal accelerations and 
eigenvalues of the tidal tensor in y-direction are equal to 
those in z-direction. In order to formulate a conservative 
criterion, we define the regions where the compression dom¬ 
inates the expansion with \a y /a x \ > 1 and A 2 /A 1 < —1. If 
the z-direction, where the tidal held is as compressive as in 
the y-direction, is taken into account as well, the above con¬ 
ditions would be \a y /a x \ > 0.5 and A 2 /A 1 < —0.5, which 
increases the distance range where the survival of TDGs is 
boosted. 

Internal kinematics: We neglect all internal kinemat¬ 
ics, such as accelerations caused by self-gravity of the stellar 
and gaseous components, kinetic stellar feedback or rota¬ 
tion. All these processes exist and their signihcance may 
vary greatly. Therefore, we do not address the very inner 
part of both TDGs and isolated DGs and focus our anal¬ 
ysis on the outer regimes where the DM or the tidal held 
dominates over the self-gravity of the baryons. Note that 
tidal compression is more important than tidal expansion, 
























Tides or dark matter sub-halos: Which ones are more attractive? 11 


as it brings material close to the centre, where it can become 
trapped by the self-gravity of the baryonic matter. For the 
detailed study of the evolution of individual TDGs or DGs, 
computationally expensive, high-resolution simulations are 
necessary. The method proposed in this work allows a quick 
way to gauge the importance of a tidal field. 

General remark: During the very early stages of the 
galaxy interaction, the rapidly changing total gravitational 
potential leads to regions of fully compressive tidal modes 
that are only short lived (Renaud et al. 20091. These re¬ 
gions can enhance the formation of bound structures. For 
later stages of the interaction, when the TDG is orbiting 
around the host galaxy, the method presented here can help 
to understand the role of the tidal field in the long-term 
survival. The two approaches are therefore complementary. 
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